Stable and Metastable Structures of Cobalt on Cu(OOl): An ab initio Study 
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We report results of density-functional theory calculations on the structural, magnetic, and elec- 
tronic properties of (1 x l)-structures of Co on Cu(OOl) for coverages up to two monolayers. In 
particular we discuss the tendency towards phase separation in Co islands and the possibility of seg- 
regation of Cu on top of the Co-film. A sandwich structure consisting of a bilayer Co-film covered 
by 1ML of Cu is found to be the lowest-energy configuration. We also discuss a bilayer c(2 x 2)-alloy 
which may form due to kinetic reasons, or be stabilized at strained surface regions. Furthermore, we 
study the influence of magnetism on the various structures and, e.g., find that Co adlayers induce a 
weak spin-density wave in the copper substrate. 
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I. INTRODUCTION 

Heteroepitaxial structures of Co and Cu ex- 
hibit intriguing magnetic properties such as n giant 
magnetoresistancea, interlayer exchange couplinga, and 
surface magnetic anisotropyu. Since these properties are 
closely related to the surface and interface morphology, 
identification and understanding of the atomic structures 
and energetics of the adsorption of cobalt on the copper 
surface are of great interest. Specifically we discuss in 
this paper the [001]-surface orientation. Thin films de- 
posited on a substrate of a different material are gen- 
erally subject to strain arising from the different lattice 
parameters of the adsorbate and substrate. Our calcu- 
lations show that the lattice constant of a ferromagnetic 
fee bulk phase of Co is 2.8% smaller than that of a fee 
Cu crystal, while the lattice constant of a hypothetical 
nonmagnetic fee cobalt crystal is 4.3% smaller than that 
of the copper crystal. Here we take the fee structure of 
cobalt, because it has been shown that a thick epitaxial 
cobalt film on Cu(OOl) can be characterized in terms of 
a tetragonally distorted face centered cubic (fct) phasea. 
The lattice mismatch between cobalt and copper sug- 
gests a small tensile strain. However for ultrathin films 
(O < 2 ML) the comparison of the bulk phases of adsor- 
bate and substrate is not necessarily very relevant. For 
example total-energy calculationsO show that the equi- 
librium lattice constant of an unsupported Co monolayer 
is 14.1% (nonmagnetic case) and 12.2% (ferromagnetic 
case) smaller than the Cu bulk lattice constant, imply- 
ing that an ultrathin film might be subject to a much 
stronger tensile strain than a thick overlayer. The re- 
lation between lattice mismatch and relaxation of the 



interlayer spacing will be discussed in Section IV below. 
While experimental studies of coverages above 2ML 
show that growth proceeds in an almost perfect layer- 
by-layer mode, for the initial two layers a deviation 
from the Frank-van der Merwe (FM) growth mode 
and a stipng dependence on the growth conditions was 
reported! Angle-resorwed X-ray photoemission spec- 
troscopy (ARXPS) dataQ indicate that the second layer 
begins to form before the first layer is completed. As- 



suming the coexistence of areas of clean Cu(OOl) surface, 
and of monolayer and bilayer islands, a LEED analysisu 
estimated that for a total coverage of one monolayer and 
deposition rates ranging between 0.016 and 0.33 ML/s 
the area covered by bilayer islands j-at room tempera- 
ture is 20 - 40%. Fassbender et alB performed STM- 
experiments at room temperature for a total coverage of 
1.35 ML and report that the fractional layer filling de- 
pends strongly on the deposition rate. For a low (0.003 
ML/s) deposition rate they found that the first layer was 
closed and 0.35ML were in the second layer, while for a 
high deposition rate (0.3 ML/s) 15% of the surface was 
still uncovered and about 50% of the surface was already 
covered by bilayer high islands. a 

X-ray photoemission scattering (XPS)LJ, AiJger elec- 
tron scattering(AES) and STM measurementsa show an 
increase of the Cu signal and decrease of the Co signal 
upon annealing which was interpreted as segregation of 
substrate material on top of the cobalt layer. Similar 
results were reported for Fe/Cu(001)t3. This effect was 
explained in terms of the lower surface energy of Cu com- 
pared to Co. We note that the application of this argu- 
ment to thin film systems is not trivial because of the 
energy cost of the additionally created Cu/Co interface. 
Yet our studies show that in the case of Co on Cu(001) 
the contribution of the interface energy is very small (sec 
Section |ffl|). 

The impact of morphological changes on the magnetic 
properties of Co/Cu(001) was recently investigated with 
X-ray magnetic circular dichpaism (XMCD) and magneto 
optical Kerr effect (MOKE)U experiments. At a cobalt 
coverage of 1.8 ML a sudden jump of the Curie temper- 
ature was measured which changed strongly with time 
or a subsequent heat treatment. The authors speculated 
that the critical thickness coincides with the thickness at 
which bilayer cobalt islands coalesce. 

Depending on growth conditions (temperature, deposi- 
tion rate), significantly different structures are observed 
experimentally. Although the magnetic properties of Co 
on Cu(001) have been the subject of many theoretical 
studies, a systematic theoretical analysis of the different 
configurations and their relative stability is still lacking. 



Moreover most of the calculationstSllI 4 ] have used slabs 
with atomic positions frozen to the bulk coordinates of 
the substrate, neglecting thus the structural relaxation 
of the clean Cu(OOl) surface and of the Co/Cu(001) ad- 
sorbate system. 

In this paper we focus on the behavior of Co on 
Cu(OOl) under thermodynamic equilibrium conditions. 
We performed density-functional theory calculations con- 
sidering a variety of configurations (9 < 2 ML). In partic- 
ular we discuss two aspects: the formation of multilayer 
cobalt islands and sandwich structures with a copper cap- 
ping. For each system we performed a full structure opti- 
mization and establish the relation between the energetic 
trends and the structural, magnetic and electronic prop- 
erties. The paper is organized as follows: The details of 
the calculations are given in Section pi In Section III we 
discuss the stability of the systems against separation in 
multilayer islands and the influence of the capping layer. 
The structural (Section IV), magnetic (Section [V|) and 
electronic (Section VI) properties of mono- and bilayer 
cobalt films on Cu(001), as well as of the correspond- 
ing cop per c apped systems are investigated. Finally in 
Section |VII| we address the similarities and differences 
between Co/Cu(001) and Co/Cu(lll) refer™? to STM 
and ab initio results for the [1 1 ll-orie ntatiorill The re- 
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suits are summarized in Section VIII 



II. CALCULATIONAL DETAILS 

Our calculations are performed using density- 
functional theory (DFT) . The exchange-correlation func- 
tional is treated within the local-density approxima- 
tion (LDA)t3 and for the magnetic systems we per- 
formed spin-polarized calculations within the local spin- 
density approximation (LSDA). We also examined the 
possible importance of non-local exchange correlation ef- 
fects by employing the generalized-gradient approxima- 
tion (GGA) in— the parameterization of Perdew, Burke, 
and ErnzerhoiO. The results show that for our study 
LDA and GGA give the same structural and energetic 
trends. More details on this issue will be discussed in the 
Appendix. 

The Kohn-Sham equation was solved applying the full- 
potentiaLlkiearized augmented plane wave (FP-LAPW) 
methocOO. The surface is simulated by repeated slabs 
separated in z-direction by a vacuum region. Co is ad- 
sorbed on both sides of the substrate. The thickness of 
the vacuum region between the slabs, corresponding to 6 
Cu layers (10.65 A), is found to be sufficient to avoid in- 
teractions of the Co atoms. The interlayer distances dn 
and (I23 were optimized with a damped Newton dynam- 
ics and the relaxations Ad^/rfo and Ac^/do are given 
with respect to the interlayer spacing of a Cu crystal, do- 
Referring the Co-Cu and for Co bilayer systems even the 
Co-Co interlayer distances to the interlayer spacing of Cu 
is probably not an optimum choice, but it is well defined 
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TABLE I. Convergence tests performed within LDA for 
a 5-layer slab of Co(001) strained at the lattice constant of 
copper and interlayer distance optimized for JVk.. = 28. The 
surface energy E and work function cj> are given as a function 
of the plane wave cut-off E cut and the number of ky -points in 
the irreducible part of the Brillouin zone JVk N . 



and has been the common practice for such adsorbate 
systems. We therefore use this convention here as well. 

The lattice constant for the fee copper crystal ac u = 
3.55 A, obtained from a non-relativistic calculation, is 
1.6% smaller than the measured one (3.61 A), 0.1% of 
which reflects our neglect of zero point vibrations in the 
theory. The lateral lattice parameter of the Cu substrate 
was set to the calculated lattice constant for a fee copper 
crystal. We chose a muffin tin (MT) radius of R^J = 
2.20 bohr for the Cu atoms and a slightly smaller radius 
Rqv — 2.15 bohr for the Co atoms to prevent overlap of 
the MT spheres due to the strong relaxation found for 
some systems. 

The stability of various systems is analyzed with re- 
spect to the formation energy. Assuming that the slab 
is in thermal equilibrium with a Co and a Cu crystal, 
acting as reservoirs of Co and Cu atoms, the formation 
energy in eV per (1 x l)-unit cell is defined as: 
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where A is. the area of the surface unit cell of the consid- 
ered slabcil and the factor 2 accounts for the presence of 
two surfaces of the slab. Aq u and Nco are the number 
of Cu and Co atoms in the slab supercell and E^ k and 
_E^" lk are the energies of a Cu or a Co atom in the respec- 
tive fee bulk crystals at the theoretical equilibrium lattice 
constants. Thus for a pure Cu slab (Nco = 0) E is the 
Cu surface energy, and for a pure cobalt slab (Nqu = 0) 
with aii = aco it is the surface energy of cobalt. 

The LAPW wave functions within the muffin tins 
(MTs) were expanded in spherical-harmonics with an- 
gular momenta up to l^ ax = 10. Non-spherical contri- 
butions to the electron density and potential within the 
MTs were considered up to l^ x = 4. The cutoff for the 
Fourier-series expansion of the interstitial electron den- 
sity and potential was chosen to be G max = 12.0 bohr -1 . 
Extensive convergence tests with respect to k|| -point set 
and the energy cutoff for the basis set were performed 
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TABLE II. DFT-LDA results for the surface energy E l , 
work functions cf>, and relaxation Adi2/do of the clean 
Cu(001)-surface are given as a function of the number of slab 
layers A^ laycr , N k = 21, E cut = 15.0 [Ry]. 
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5.31 
5.30 
5.30 



TABLE III. Formation energy and work functions for 
lCo/Cu(001) with optimized interlayer distances as a func- 
tion of the number of substrate layers Mayer. 



for a 5-layer Co(001) slab at the lattice constant of cop- 
per and relaxed interlayer distance. The results are 
shown in Table |. A numerical accuracy of 6% for 
the formation energy is achieved with -Ecut = 12.8 Ry, 
while Ecut = 15.6 Ry is needed for an accuracy of 1%. 
Thus a cutoff parameter of 15.6 Ry was chosen through- 
out the calculations. The Brillouin-zone integration was 
performed with a special point .-set generated after the 
scheme of Monkhorst and Packed We obtained an ac- 
curacy of the Brillouin-zone integration better than 1% 
by using 21 kii-points in the irreducible wedge of the 
Brillouin-zone (IBZ) (see Table |). 

The bulk energies needed as a reference to determine 
the formation energy (see Eq.(Jl|)) were calculated using 
the same LAPW-parameters as in the slab calculations. 
For the bulk calculation 104 k-points in the IBZ were 
used. 

Prior to investigating the effects of adsorption of 
cobalt, we checked the required thickness of a copper 
slab, to ensure a good representation of the properties 
of the clean Cu(001)-surface. The surface energies, work 
functions and interlayer relaxations for slabs of 3, 5, 7, 
and 9 layers of copper are compared in Table |l|. For a 
5-layer slab the relaxation between surface and subsur- 
face layer of — 3.0L% is close to the experimental value 
obtained by MEISE3 (-2.4%) while the LEED-resultN 
is smaller (—1.1 ± 0.4%). The calculated work func- 
tion </> = 4.78 eY is in good agreement with experiment: 
4.59 ± 0.05 eVEa, 4.76 eVB, and 4.77 ± 0.05 eVEI. 

A further requirement of the thickness of the slab is 
that the interaction of the Co layers on both sides of the 
slab trough the substrate is negligible for the questions 
of concern. To test the strength of this interaction we 
studied the formation energy and work function for 1ML 
Co on Cu(001) (lCo/Cu(001)) as a function of the sub- 
strate thickness. The results are summarized in Table 
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TABLE IV. Formation energies E { and the work functions 
for various structures. 



III. Both the formation energy and work function con- 
verge quickly with the substrate thickness. On the basis 
of these results we conclude that a 5-layer copper slab 
represents a good approximation of the Cu(001)-surface. 



III. FORMATION ENERGY AND STABILITY 

In order to identify the equilibrium configuration of 
Co on Cu(001) we investigate in this Section the ten- 
dency towards separation in multilayer islands and the 
influence of a copper capping layer. The studied sys- 
tems include the clean Cu(001) surface, a monolayer 
and bilayer thick cobalt film on Cu(001) denoted by 
lCo/Cu(001) and 2Co/Cu(001), respectively, as well as 
the corresponding capped systems, 1Cu/1Co/Cu(001) 
and lCu/2Co/Cu(001). Additionally we investigated a 
bilayer Co-Cu-c(2 x 2)-alloy. The calculated formation 



energies and the work functions are given in Table IV . We 



consider several ways in which a total coverage of 1ML 
Co can be arranged on a Cu(001) surface. The energy of 
a system consisting of more than one domain, namely re- 
gions of clean copper surface and regions covered by large 
cobalt islands, is simply given by the weighted sum of the 
formation energy of the clean Cu(001) surface and those 
of the Co-island. Under the assumption that the islands 
are large, the contributions of the step edges and side 
facets of the islands are negligible and were not taken into 
account. A schematic presentation of the different struc- 
tures are given in Fig. [j\ together with energy changes 
with respect to the case where the whole surface is cov- 
ered by a monolayer thick (1 x l)-cobalt layer. 

Our calculations show that a monolayer film, Fig. |l](a), 
would separate into a clean Cu(001)-surface and a bilayer 
island, Fig. 0(b). For the nonmagnetic case the gain in 
energy is AE NM = 0.59 eV/(l x l)-cell and for the fer- 
romagnetic case it is AE FM = 0.41 eV/(l x l)-cell. This 
result can be explained in terms of the higher coordina- 
tion of the cobalt atoms in the bilayer film and correlates 
with the substantial broadening of the cobalt d-band and 
the strong relaxation between Co layers in 2Co/Cu(001) 
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a) close packed monolayer high 
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b) close packed bilayer high Co-island ae = 0.59 eV 
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c) bilayer Co-Cu-c(2x2)-alloy 
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AE =0.28eV 
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AE =0.15eV 



d) close packed bilayer high 
Cu+Co-island 



NM 

AE =0.49eV 

FM 

AE =0.33eV 



e) close packed trilayer high 
Cu+Co+Co-island 



NM 

AE =0.79eV 

FM 

AE =0.56eV 



FIG. 1. Schematic diagram of the different adsorbate con- 
figurations for a Co coverage of Q = 1.0 ML. The formation 
energy changes for configurations b)-e) in the nonmagnetic 
A_B NM and ferromagnetic case A_E FM is given with respect 
to formation energy of the close packed monolayer high is- 
land shown in a). 



as will be discussed later in this paper. Concerning the 
effect of magnetism, we see that it reduces, but does not 
change the tendency towards formation of bilayer islands. 
Experimental studiesQu show that copper segregates 
onto the surface after annealing. Therefore we study 
here the influence of a copper capping layer on stabil- 
ity. Covering lCo/Cu(001) with a monolayer of copper, 
Fig. [l](d), reduces the energy of the system by AE NM = 
0.49 eV/(l x l)-cell. Compared to the cobalt termi- 
nated systems the copper-capped systems gain less spin- 
polarization energy because of the hybridization with 
the capping layer. Consequently the energy gain due 
to a capping layer for the magnetically ordered system 
is lower than the one for the nonmagnetic: (AE FM — 
0.33 eV/(l x l)-cell). The influence of the capping layer 
on the magnetic properties of the copper covered sys- 
tems rWjill be discussed in Section [V| Saul and Weiss- 
mani£3 recently calculated the surface segregation en- 
ergy of 3d-impurities (Fe, Co, Ni) in Pd, Ag, and Cu. 
They found, in agreement with our results, that embed- 
ding in the bulk of the host material is connected with 
a substantial gain in energy both for nonmagnetic and 
magnetic impurities/ (layers), the effect being weaker for 
the latter. We also note that the copper capping layer 
in the 1Cu/1Co/Cu(001) and ICu/2Co/Cu(001) systems 
has properties similar to the clean Cu(001)-surface; for 
example we find that the work functions of the systems 
are (fcu(ooi) = 4.78 eV, 0icu/ico/cu(ooi) = 4.89 eV, and 



->lCu/2Co/Cu(001) 



4.74 eV. 



Recent combined STM and RHEED cxperimentsLJ de- 
tected ordered c(2 x 2) regions when a total coverage of 
1ML Co was deposited on Cu(001) at room tempera- 
ture and subsequently annealed at 450 K. Motivated by 
these results, we studied a configuration, where starting 
from lCo/Cu(001) every other Co atom is replaced by 
a Cu atom in the substrateJayer underneath. In this 
way a bilayer c(2 x 2)-alloyE3, shown schematically in 
Fig. |l](c), is formed. We find that this configuration is by 
0.28 eV/(l x l)-cell (nonmagnetic) and 0.15 eV/(l x 1)- 
cell (ferromagnetic case) more favorable than the (1 x 1)- 
monolayer in Fig. 13(a) . However, it is a metastable 
structure because transition into a cobalt monolayer cov- 
ered by copper, 1Cu/1Co/Cu(001) in Fig. 0(d), leads to 
an energy gain of 0.22 [eV/(l x l)-cell] (nonmagnetic) 
and 0.18 eV/(l x l)-cell (ferromagnetic case). Thus, 
the bilayer c(2 x 2)-alloy lies energetically between the 
lCo/Cu(001) and 1Cu/1Co/Cu(001) systems and it may 
be stabilized kinetically. This surface alloy might also 
represent a favorable configuration with respect to sur- 
face strain relief. Indeed a c(2 x 2)-pattern was observed 
experimentally preferentially in the middle of large is- 
landsEl. 

We also studied whether the cobalt (1 x l)-layer will 
prefer to be buried deeper in the substrate. Our calcula- 
tions for the 1Cu/1Co/Cu(001) and 2Cu/lCo/Cu(001) 
systems show that there is no additional energy gain 
through covering the system with a thicker copper layer. 



System 
. Cu/lCo/Cu. 
. Cu/2Co/Cu. 
. Cu/3Co/Cu. 



Eisn 



E-pt 



0.163 


0.113 


0.044 


0.066 


-0.005 


0.046 



TABLE V. Interface energies given in [eV/(l x l)-cell] for 
nonmagnetic (Elj M )a,nd ferromagnetic (-E"pm) systems as a 
function of the cobalt interlayer thickness. 



The segregation of Cu on the surface is typically ex- 
plained by the lower surface energy of Cu(001) com- 
pared to Co(001). Still this argument is only appli- 
cable if the interface energy were small and thus neg- 
ligible. In order to calculate the energy cost to cre- 
ate an interface we studied three different systems with 
one, two and three cobalt interlayers in copper bulk, 
marked as (. . . Cu/lCo/Cu. ..)>(••• Cu/2Co/Cu. . . ) and 
(. . . Cu/3Co/Cu. . . ), respectively. Here the lateral pa- 
rameter is set to the copper bulk value while the in- 
terlayer distances are relaxed. The interface energies, 
E l , are calculated analogously to the formations energies 
[compare Eq.Q]. In order to subtract the effect of elastic 
strain due to the lattice mismatch of the two materials 
we use as a reference energy for cobalt the bulk energy 
of fct cobalt with an = acu and relaxed a z c3 instead 
of the one of fee cobalt at the cobalt lattice constant. 
Table |V| lists the results for nonmagnetic and ferromag- 
netic cobalt interlayers. It shows that the values of E 1 
for a cobalt bilayer and trilayer are very close, i.e. the 
interface energy converges quickly with the thickness of 
the cobalt layer. The interface energy is indeed signif- 
icantly smaller than the difference of the surface ener- 
gies of the clean Cu(001)-surface {E f = 0.79 eV/(l x 1)- 
cell) and a thick Co(001) film with a lateral parameter 
fixed to the lattice constant of copper and relaxed in- 
terlayer distances (E ! NM =7-4 -21 eV/(l x l)-cellEl and 
Ep M = 1.11 eV/(l x l)-cellE3). For this reason the com- 
mon argument that simply the surface energy difference 
of cobalt and copper explains the segregation of sub- 
strate material on the surface works in the case of Co 
on Cu(001). 

In analogy to the Co terminated system, the single 
Co layer capped by Cu shown in Fig. [|(d) will tend 
to separate into a clean Cu(001) surface and a double 
Co layer capped by Cu, Fig. 0(e). Still the energy gain 
due to phase separation [AE™ = 0.31 eV/(l x l)-cell, 
AE FM = 0.23 eV/(l x l)-cell)] is only about half the en- 
ergy gain for the system with Co on the surface. We can 
summarize that both the magnetic ordering and the cap- 
ping layer weaken the tendency towards cobalt clustering 
but qualitatively we observe the same behavior with and 
without magnetism. 

Yet we need to find out whether the bilayer film will 
be stable or if higher cobalt islands may form. In order 
to determine the formation energy of a ./V-layer (TV > 2) 
thick cobalt island we assume that each of the interme- 
diate cobalt layers has an energy of a bulk atom in a 
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) uncapped Co-islands 
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2 4 6 8 

Thickness of Co-islands [ML] 

FIG. 2. Formation energy of different ferromagnetically 
ordered configurations for a total cobalt coverage of 1 ML as 
a function of the cobalt island thickness N. The structures 
consist of clean Cu(001) and a compact island with N Co 
layers (o) or N Co-layers capped by copper. The area covered 
by the cobalt islands is i of the whole surface. Especially for 
the copper terminated systems the separation in higher than 
bilayer cobalt islands is unlikely because of a negligible energy 
gain. 



tetragonal cobalt crystal with a» — acu and fully relaxed 
a 2 E3. The elastic energy contribution is the difference 
between the energy of cobalt bulk at the fee cobalt lat- 
tice constant and the energy of fct cobalt as described 
above and amounts to 0.11 eV (NM) and 0.08 eV (FM) 
per cobalt atom. Thus the formation energy of a JV-layer 
cobalt film is 
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last. 



(2) 



For a total coverage of 1ML the formation energy of a 
configuration consisting of clean Cu(001) surface and a 
TV-layers high cobalt island is given by: 
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An analogous expression for the formation energy holds 
for the copper capped systems. The formation energy 
of the ferromagnetically ordered capped and uncapped 
systems is plotted in Fig. g as a function of cobalt-island 
height N. In the following we will concentrate on the Cu- 
terminated systems because, as can be seen from Fig. ||, 
they are always lower in energy. The substantial energy 
gain due to separation of a monolayer-thick cobalt film in 
bilayer islands was already discussed above. Yet further 
separation in higher cobalt islands brings only a small 
energy gain, e.g., the gain due to separation from bilayer 
in trilayer islands for the Cu-terminated islands is about 
0.03 eV/(l x l)-cell. We note that this energy gain is 
mainly due to the increase of the clean Cu(001)-surface. 
Moreover the cost of the island facets, which was not 
taken into account in the present discussion, grows with 
island height. Because of the small energy gain, as shown 
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FP-LAPW NM 
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-7.0% 
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-3.0% 
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1ML Cu/2ML Co/Cu(001) 


FP-LAPW NM 
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-5.0% 
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-14.8% 
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TABLE VI. Relaxation Adi2/do and Adn/do of the inter- 
layer spacing in % for the first two layers compared to the 
lattice parameter of Cu bulk, do. 



in Fig. g, and the increasing cost of the sidewalk the for- 
mation of islands higher than bilayer is unlikely. We con- 
clude that the ferromagnetically ordered configuration in 
Fig. 0(e), which is by 0.56 eV/(l x l)-cell more favorable 
than the one in Fig. |l](a), represents the thermodynamic 
equilibrium structure. 



IV. STRUCTURAL PROPERTIES 

In the previous Section we identified the bilayer cobalt 
island covered by a copper capping layer as the thermo- 
dynamically stable structure. However, crystal growth 
represents a situation which is more or less far from 
thermodynamic equilibrium therefore not only the equi- 
librium structure but also other, metastable, structures 
(e.g. those shown in Fig. \j\) may occur. In this Section we 
present the results of a geometry optimization for mono- 
layer and bilayer (1 x 1)- as well as copper-capped sys- 
tems. The Co films are assumed to grow pseudomorphi- 
cally on the Cu(001)-surface, adopting the lateral spac- 
ing of the Cu crystal. The calculated relaxation of the 
interlayer spacing for the nonmagnetic and ferromagnetic 
systems is given in Table VI. We remind the reader that 
all relaxations are given with respect to the interlayer 
spacing in copper bulk. 

The first interlayer spacing in the monolayer film, 
rfco-Cu, shows an inward relaxation of 4.7% for the non- 
magnetic case, which reduces to 3.0% for the ferromag- 
netic film. At the same time the interlayer spacing be- 
tween the first and second substrate layer, which is con- 
tracted by 3.0% for the clean Cu(001) surface, expands 
back to the bulk value —0.3% (0.0%) for the nonmagnetic 
(ferromagnetic) system. This result can be explained 
in terms of the bond-cutting model. Due to the miss- 



ing bonds of the surface atoms the strength of the re- 
maining bonds to the subsurface-Layer is enhanced, giv- 
ing rise to an inward relaxationEl The.-bond strength 
is also related to the d-band occupationcHl, thus Co-Cu 
bonds are stronger than Cu-Cu bonds and consequently 
upon cobalt adsorption we observe a stronger relaxation 
of the Co-Cu-interlayer distance, while the Cu-Cu dis- 
tance expands. Previous ab initio resultsEj found a re- 
laxation of the first interlayer distance of the ferromag- 
netic lCo/Cu(001) surface of -10.4%. The reason for 
the discrepancy with our result (which is —3.0%) is in the 
choice of the lateral lattice parameter. While we use a 
non-relativistic treatment of the valence electrons and the 
corresponding theoretical equilibrium lattice constant of 
Cu (3.55 A), Wu and FreemanES used a semi-relativistic 
treatment which gives a noticeably smaller lattice con- 
stant (3.52 A). Nevertheless, in their adsorbate study 
Wu and Freemanll3 set the lateral lattice parameter to 
the substantially larger experimental value (3.61 A). As 
a consequence, the, strong interlayer relaxation found by 
Wu and FreemanEj just reflects that their copper surface 
is under tensile strain. We tested this and indeed could 
reproduce the effect: When we use a semi-relativistic 
treatment of the valence electrons and still force the Cu 
substrate to assume the experimental lattice constant we 
obtain Ad 12 /d = -7.9% and Ad 23 /d = -2.3%. 

For the bilayer Co film we obtain a surprisingly strong 
contraction of the interlayer distance of c?co-Co of —17% 
in the nonmagnetic case. For the ferromagnetically 
ordered system the contraction is somewhat smaller, 
— 13.4%, due to the magnetovolume effect. These re- 
sults can hardly be explained by comparing the lattice 
constant of the fee bulk phase of Co with that of bulk 
Cu. Such comparison would give a lattice mismatch of 
—4.3% in the nonmagnetic and —2.8% in the ferromag- 
netic case. Thus, in such description one would say that 
the Co film is strained, but that the effect is not very 
large. However, it is questionable whether the compari- 
son of the bulk lattice parameters of the two materials, is 
a good approach for understanding ultrathin films with 
O < 2 ML. For example, tatal-energy calculations based 
on the FP-LAPW methodu show that the difference be- 
tween the equilibrium lattice constant of a free standing 
cobalt monolayer and that of the Cu substrate is —14.1% 
for a-annmagnetic and —12.2% for a ferromagnetic mono- 
layerEl Therefore, if we refer the strain in the Co adlay- 
ers to the lattice parameter of the free-standing Co layer, 
the strain is significant, and the above noted interlayer 
relaxation then simply reflects the reaction of the Co film 
to this big strain. Indeed, we think that this description 
is appropriate (in a qualitative sense) because for a very 
thin cobalt film (9 < 2 ML) the bonding to the noble 
metal substrate can only partially replace the bonds to 
missing cobalt neighbors. Thus, the adsorbed film will 
still bear some resemblance to the free-standing one. The 
above result also indicates that the weaker binding to the 
substrate is balanced by forming a strong bond between 
the two cobalt layers. 



The competition between Co-Co and Co-Cu-bonding 
is also a driving force for the structural changes in the 
capped systems. The hybridization with the copper cap- 
ping layer has the general effect of weakening the exist- 
ing Co-Cu and Co-Co bonds in the lCo/Cu(001) and 
the 2Co/Cu(001) system, respectively. Consequently, 
the interlayer distance between the cobalt and copper 
layer increases from -4.7% (NM) and -3.0% (FM) 
in lCo/Cu(001) to -3.0% (NM) and -2.0% (FM) in 
1Cu/1Co/Cu(001). Similarly the strong relaxation be- 
tween the two cobalt layers decreases from —17.0% (NM) 
and -13.4% (FM) in 2Co/Cu(001) to -14.8% (NM) and 
-12.5% (FM) in lCu/2Co/Cu(001). On the other hand 
the stronger Co-Co bond induces a weaker binding with 
the capping layer which is reflected in the smaller inward 
relaxation of the distance between the capping layer and 
the cobalt film of -5.0 % for the lCu/2Co/Cu(001) sys- 
tem compared to -7.0% for the 1Cu/1Co/Cu(001) sys- 
tem. 



Table VI also contains structural data determined with 



LEED. We note, however, that such structural analysis 
is complicated and not unambiguous, because, as dis- 
cussed above, in the Co/Cu(001) system several domains 
and/or metastable structures may coexist. In the ab- 
sence of knowledge about these various structures and 
their energies^ it appeared to be a reasonable choice for 
Clarke et alE3 to assume that Co on Cu(001) will form 
a full Co monolayer. And based on this assumption 
they determined an inward relaxation of —6% for both 
dco-Cu an d rfcu-Cu- Our work, however, shows that the 
lCo/Cu(001) system is unstable with respect to the for- 
mation of bilayer islands and capped bilaygrs structures. 
In a more recent LEED study Cerda et alB assumed the 
coexistence of regions of clean Cu(001), Co monolayer 
and Co bilayer islands, but Co layers with a Cu capping 
layer, which we find to have the lowest total energy, were 
not considered. Thus, in both experimental analyses the 
model assumptions didn't include all relevant systems. 



V. MAGNETIC PROPERTIES 

The layer-resolved magnetic moments in the four sys- 
tems studied, lCo/Cu(001), 
2Co/Cu(001), 1Cu/1Co/Cu(001) and lCu/2Co/Cu(001) 
are given in Table VII. To be precise, these are the 
contributions from the muffin tin region only. The top 



layer in lCo/Cu(001) exhibits an enhanced magnetic mo- 
ment (M Co ( S ) = 1.71 /is) compared to the bulk value of 
1.52 /is, calculated at the equilibrium lattice constant of 
cobalt. This is due to the larger lateral constant of the 
epitaxial cobalt adlayer and to the reduced coordination 
on the surface. Further we find that the surface layer 
of the 2Co/Cu(001) system exhibits the same magnetic 
moment (1.71 /is) as the lCo/Cu(001) system. In fact, 
a thick fee cobalt film at the lattice constant of copper 
also has a similar moment, namely 1.78 /is. However, 



System 


layer 


M [fi B ] 


lCo/Cu(001) 


Cu(C) 
Cu(S-2) 
Cu(S-l) 

Co(S) 


-0.004 
-0.014 
0.024 
1.711 


2Co/Cu(001) 


Cu(C) 
Cu(S-3) 
Cu(S-2) 
Co(S-l) 

Co(S) 


-0.002 
-0.009 
0.016 
1.472 
1.706 


Co(001) at ac„ 


Co(C) 

Co(S-l) 

Co(S) 


1.648 
1.615 
1.783 


1Cu/1Co/Cu(001) 


Cu(C) 
Cu(S-3) 
Cu(S-2) 
Co(S-l) 

Cu(S) 


-0.001 
-0.008 
0.027 
1.445 
0.040 


lCu/2Co/Cu(001) 


Cu(C) 
Cu(S-4) 
Cu(S-3) 
Co(S-2) 
Co(S-l) 

Cu(S) 


-0.001 
-0.007 
0.022 
1.383 
1.374 
0.035 



TABLE VII. Layer resolved local magnetic moments in the 
ferromagnetic systems in [/is] as obtained from the slab cal- 
culation. S, S-l, S-2, etc. denote the position of the corre- 
sponding layer with respect to the surface, S being the surface 
layer and C being the central layer of the slab. 



the magnetic moment of the subsurface Co layer, which 
binds to the Cu substrate, is reduced to 1.47 /is- The 
corresponding magnetic moment of subsurface cobalt in 
a thick fee cobalt film at the lattice constant of copper 
is 1.62 /is. The lower magnetic moment of subsurface 
cobalt is a consequence of the higher coordination and 
the strong contraction of the interlayer spacing dco-Co- 

The hybridization with the copper capping layer re- 
duces the magnetic moment of the first Co layer both 
in the 1Cu/1Co/Cu(001) and lCu/2Co/Cu(001) sys- 
tems by about 0.3 /is compared to the lCo/Cu(001) and 
2Co/Cu(001) systems. It is interesting to note that both 
Co layers in lCu/2Co/Cu(001) have the same magnetic 
moment (1.38 /is) which can be explained by the fact 
that Co(S-l) and Co(S-2) have the same coordination of 
Co and Cu atoms. 

To our knowledge magnetic moments for lCo/Cu(001) 
have not yet been measured due to the already dis- 
cussed difficulties in the preparation of a single cobalt 
monolayer on Cu(001). Our calculated value of the 
surface magnetic moment of lCo/Cu(001), 1.71 /is, is 
slightly lower than, that obtained, in previous calcula- 
tions, e^ 1.78/isli3 and 1.76 hJB from FP-LAPW and 
1.85/islij from FP-LMTO (full-potential linearized muf- 
fin tin orbitals) calculations. The differences are at- 
tributed mainly to the use of the experimental lattice 
constant of copper (3.61 A) and/or the lack of consider- 



ing the inter layer relaxation in Refs.li3 EJ (see also Sec- 
tion 0). 

The magnetic moments of 1.9 ML and 2.1 ML of 
Co deposited on Cu(OOl) measured with x-ray-jpag- 
netic circular jdichroism (XMCD) are 1.71 ±0.1 /isEII and 
1.77 ± 0.1 /i B H respectively. The XMCD-spectra were 
recorded at 40 K, but information about the prepara- 
tion conditions, which could tell whether the Co layers 
were capped by Cu, is not available. Yet, the mag- 
netic moment compares well with our calculated mag- 
netic moment for 2Co/Cu(001). With respect to other 
theoretical work, we note that the same trend of an en- 
hanced magnetic moment in the surface layer (1.85 /is) 
and a reduced magnetic moment in the subsurface layer 
(1.60 hb) was fou-nd in a previous FP-LAPW calculation 
for 2Co/Cu(001)c3. In this study the lateral parameter 
was fixed to the experimental lattice constant of copper 
and relaxation of the interlayer spacing was not taken 
into account. However, the strong relaxation of the in- 
terlayer spacing in 2Co/Cu(001) discussed in Section IV 
has a noticeable influence on the magnetic moments and 
cannot be neglected. 

Our results reveal also that the adsorbed cobalt film in- 
duces a small polarization in the substrate. The magnetic 
moment of the copper layer at the interface is positive, 
e.g. in lCo/Cu(001) it is 0.024 \i B . Then, in the next 
layer it switches to a negative value (— 0.014 \ib). Also 
the central layer of our 5-layer Cu slab has a very small 
negative moment, —0.004 jib- The oscillation of the mag- 
netic moment perpendicular to the surface indicates the 
formation of a spin-density wave. This striking effect is 
observed for all four studied systems. However, we note 
that for a detailed investigation of this effect, a thicker 
substrate slab has to be considered. The magnetic mo- 
ment induced in the capping layer is somewhat larger 
than the one induced in the substrate layer: 0.040 Hb m 
lCu/lCo/Cu(001) and 0.035 \i B in lCu/2Co/Cu(001). 



VI. ELECTRONIC PROPERTIES 

The calculated electronic properties are consistent with 
the above discussed structural and energetic trends. Fig- 
ure p| shows the local density of states (LDOS) of the 
d-bands of the adsorbate and substrate layers obtained 
from a nonmagnetic calculation. For lCo/Cu(001) the 
Co d-band is rather narrow, the LDOS at the Fermi level 
is very high and the overlap with the copper d-band is 
small, which reflects that the interaction between Co and 
Cu is not very strong. For the 2Co/Cu(001) system the 
d-states of the surface and subsurface Co layers overlap 
and their d-bands receive a substantial broadening. At 
the same time the LDOS at the Fermi level is lowered. 
The broadening of the cobalt d-bands in 2Co/Cu(001) is 
an indication for the strong interaction between the two 
cobalt layers. The same effect of broadening of the ri- 
band of Co is observed for lCu/2Co/Cu(001) compared 



to the Co d-band in 1Cu/1Co/Cu(001). 

The layer-resolved LDOS of the d-bands for the fer- 
romagnetic systems is given in Fig. 0. The major- 
ity band of Co is completely filled and the minor- 
ity band is only partly filled, reflecting the fact that 
Co is a strong ferromagnet. For 1CoFM/Cu(001) and 
1Cu/1CoFM/Cu(001) the Fermi level crosses the minor- 
ity d-band of cobalt almost at its maximum while for 
2CoFM/Cu(001) and 1Cu/2CoFM/Cu(001) the Fermi 
level lies in a dip of the Co d-bands. A "harder" elec- 
tronic structure, i.e. lower density at the Fermi level is 
typically considered an indication for a more stable sys- 
tem. 

Both majority and minority d-bands of copper are oc- 
cupied and lie ca. 2 eV below the Fermi level. Still the 
minority and majority d-band have a very different struc- 
ture with the majority band being broader in general. 
Actually the band width correlates with the strength of 
interaction with the cobalt film: While there is a substan- 
tial overlap between the majority d-bands of cobalt and 
the substrate layer beneath or the capping layer above, 
the corresponding minority bands have a very small over- 
lap. 



VII. COMPARISON OF CO/CU(001) AND 

CO/CU(lll) 

Prior to our work Pedersen et al&3 studied the growth 
of Co on the (111) surface of Cu with STM and LMTO- 
calculations. The STM measurements showed that the 
islands consist of several cobalt layers with the lowest 
layer possibly growing subsurface. At elevated temper- 
atures vacancy islands formed in the terraces close to 
steps and the substrate material etched from these holes 
covered the cobalt islands. 

Comparing the LMTO calcinations for nonmagnetic 
systems with [lll]-orientatior£3 with our FP-LAPW re- 
sults for the [001]-orientation we note that the general 
behavior is similar: The systems are unstable against 
phase separation and clustering. The energy gain from 
the separation of a monolayer film in a bilayer island 
and a clean Cu surface is A-Emoi) = 0.59 eV/(l x 1)- 
cell and A£/m) = 0.39 eV/(l x l)-cell respectively. 
The corresponding energy gain for the capped systems 



is AE, 



(100) 



0.31 eV/(l x l)-cell and AE (111 



(in) 



0.18 eV/(l x l)-cell. 

The segregation of Cu onto the surface lowers the 
energy of the system for both orientations: for mono- 
layer coverages AE^qo) = 0.49 eV/(l x l)-cell and 
AE (111) = 0.30 eV/(l x l)-cell; for 2ML of Co AE {100) = 
0.42 eV/(l x l)-cell and AE {111) = 0.20 eV/(l x l)-cell. 

Still in all cases the energy gain is lower for the (111) 
surface than the (100) surface. This trend reflects the 
difference in coordination numbers: In a bond cutting 
model of metallic bonding the energy of an atom roi. 
scales as the square root of the local coordinationol 
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FIG. 3. Local density of states of the d-bands of the dif- 
ferent atomic layers in the nonmagnetic systems. S, S-l, S-2, 
S-3 denote the surface and the subsequent subsurface layers 
respectively. We display the contribution from inside the muf- 
fin-tin spheres. The Co bands are represented with dashed 
and long-dashed lines, the Cu bands with a solid line. The 
calculated LDOS was broadened by a Gaussian with a width 
of 2a = 0.2 eV. 
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FIG. 4. LDOS of the d-bands of the different atomic layers 
in the ferromagnetic systems. S, S-l, S-2, S-3 denote the sur- 
face and the subsequent subsurface layers respectively. The 
contribution from inside the muffin-tin spheres is displayed. 
The Co bands are marked with dashed and long-dashed lines, 
the Cu bands with a solid line. The calculated LDOS was 
broadened by a Gaussian with a width of 2<r = 0.2 eV. 



Adsorption of a Co layer or of a Cu capping layer implies 
a change of coordination of the atoms in the added layer 
from 4 to 8 for the [100]- and from 6 to 9 for the [111]- 
orientation. And for the atoms in the layer, which after 
adsorption becomes the second layer, the coordination 
changes from 8 to 12 for the [100]- and from 9 to 12 for 
the [lll]-orientation. Thus, the energy gain is smaller 
for the [lll]-orientation than for the [100]-orientation. 



VIII. SUMMARY 

In summary we identify a bilayer cobalt island covered 
by copper as the lowest energy configuration. However, 
growth is ruled by kinetics. Therefore, it is to be expected 
that under realistic conditions metastable structures may 
exist at surfaces, and some examples were identified in 
this paper. Total-energy considerations show that the 
(1 x l)-film tends to separate into areas of bilayer cobalt 
islands and clean copper surface, and this is indeed \l 
with experimental observations of bilayer growtl 
Our total energy and electron density of states results 
show that the stability of the bilayer film is due to the fact 
that the Co atoms prefer to attain a high coordination 
of alike atoms. A consequence of this is the very strong 
contraction of the interlayer distance between the two 
cobalt layers and the substantial broadening of the cobalt 
d-bands in the adsorbed Co-bilayer as compared to the 
d-band of a single Co adlayer. 

The segregation of substrate material onto the Co ad- 
layer results in a substantial energy gain 0.5 eV/(l x 1)- 
cell). We also studied a two-layer surface alloy of Co and 
Cu with a c(2 x 2) periodicity. This is found to be energet- 
ically less favorable than a separation into Cu-capped Co 
bilayer-adsorbates, but at strained regions of the surface 
this surface alloy may be stabilized. Indeed, a c(2 x 2) 
surface structure was observed in recent combined STM 
and RHEED experimentseia. 

Generally, the ferromagnetically ordered systems are 
lower in energy than the nonmagnetic, but the relative 
stability of different configurations remains qualitatively 
unchanged by magnetism and the structural trends are 
well described by the nonmagnetic systems. 

For low coverages (0 < 0.25 MLj , w e also find that 
cobalt may adsorb substitutionallyL-ilQ. However, with 
increasing coverage the substitutional adsorption be- 
comes energetically unfavorable compared to the forma- 
tion of compact islands. 
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System 


-E-LDA 


^GGA 


Cu(001) 


0.78 


0.61 


lCo/Cu(001) NM 
lCo/Cu(001) FM 


1.75 
1.54 


1.47 
1.22 


2Co/Cu(001 NM 
2Co/Cu(001 FM 


1.55 
1.48 


1.27 
1.11 


1Cu/1Co/Cu(001) NM 


1.26 


0.99 


lCu/2Co/Cu(001) NM 


1.13 


0.85 



TABLE VIII. The formation energies i?£ DA and Eq GA of 
the different configurations calculated within LDA and GGA, 
respectively, given in [eV/(l x l)-cell]. The lateral parameter 
is set to the corresponding (LDA or GGA) equilibrium lattice 
constant of copper. 



IX. APPENDIX 

For ses^eral systems we performed calculations with 
the LDALL3 and with the GGAEZI. The formation ener- 
gies for lCo/Cu(001) (NM and FM), 2Co/Cu(001) (NM 
and FM), 1Cu/1Co/Cu(001), and lCu/2Co/Cu(001) are 
given in Tabic VIII . The lateral parameter was set to the 
lattice constants of copper obtained within the LDA and 
GGA approach, respectively. The LDA value 3.55 A is 
1.7% smaller than the measured one, 3.61 A, while with 
the GGA the lattice parameter (3.65 A) is 1.1% bigger 
than the experimental value (zero-point vibrations are 
neglected in the theory). 

The formation energies obtained with the GGA are 
generally lower than the LDA results and the differences 
are between 0.2 and 0.3 eV/(lx)-cell. This effect was 
also observed previously for the clean copper surfaced. 
Yet the trends between the different configurations re- 
main unchanged. For example the energy gain from the 
separation of a monolayer cobalt film in a bilayer cobalt 
island and clean Cu(001)-surface is 0.60 eV/(l x l)-cell 
(LDA) and 0.53 eV/(l x l)-cell (GGA) for the non- 
magnetic systems and 0.41 eV/(l x l)-ccll (LDA) and 
0.36 eV/(l x l)-cell (GGA) for the ferromagnetically or- 
dered systems. The equilibrium configuration of clean 
copper surface with bilayer cobalt islands, covered by 
copper (see Fig. ffl(f)), is by 0.79 eV/(l x l)-cell (LDA) 
and by 0.74 eV/(l x l)-cell (GGA) more favorable than 
lCo/Cu(001) in Fig. @(a). 

A structural opti mizati on was performed for all sys- 
tems listed in Table |VII]| using both approaches, LDA 
and GGA. No noticeable differences were obtained ex- 
cept for the systems, containing a bilayer cobalt film, 
where the contraction of the distance between the 
two cobalt layers was slightly stronger with the GGA, 
e.g. for 2Co/Cu(001) Ad|£ A Co /4 DA = 17% and 
Ad^^ Co /d^ GA = 18.6 %. However, these minor differ- 



ences do not alter the discussion in Section IV 



The larger lateral parameter in GGA produces a sub- 
stantial enhancement of the magnetic moments, e.g. 
the surface magnetic moment of cobalt in lCo/Cu(001) 
changes from 1.71 /j,b (LSDA) to 1.86, /is (GGA) and in 



10 



2Co/Cu(001) from M^da = ln /i£ , and M^f s A 1} = 
1.47 M b to MgGA = L81/lB and AfGGA^ = l64mb . 
This result is not surprising and is in line with the 
changes of the magnetic moment for fee cobalt bulk from 
1.52 mb (LSDA) to 1.69 ms (GGA). 

In conclusion, both approximations of the exchange- 
correlation potential, LDA and GGA, lead to the same 
results for the structural, energetic, and magnetic prop- 
erties of the configurations studied in this work. 
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